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ABSTRACT 

An  efficiency  analysis  of  all  known  algorithms  for  finding  all 
circuits  of  a  graph  is  made,  and  bounds  on  the  computation  time  are  derived. 
Best  algorithms  are  exhibited.   The  vector-space  method  of  circuit  finding  is 
discussed  at  length.   Proven  is  the  fact  that  the  graphs  K  ,  K.  ,  K.  -x  and 

J  4      4- 

K    are  the  only  ones  which  have  2^-1  circuits  (ju  is  the  nullity  of  the 
graph),  or,  equivalently,  have  no  two  edge-disjoint  circuits.   A  class  of 
graphs,  circuit -graphs,  derivable  out  of  a  graph  are  defined.   If  an 
efficient  algorithm  for  the  enumeration  of  all  connected,  induced  subgraphs 
is  known,  these  circuit -graphs  can  advantageously  be  used  in  the  generation 
of  circuits. 
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1.   Introduction 

In  many  applications  of  graph  theory,  it  is  desired  to  obtain  a 
list  of  all  circuits  of  a  graph.   Examples  are  program  optimization,  program 
segmentation,  testing  graphs  for  isomorphism,  signal  flow-graph  theory,  etc. 
(See  Prabhaker  (1972)).   We  will  not,  however,  dwell  on  the  applications  here. 
The  problem  considered  is:  Given  an  undirected  graph  (for  definitions  of 
the  terms,  see  below),  find  all  its  circuits.   The  analogous  problem  for  a 
digraph  is  to  find  all  its  dicircuits.   The  two  problems  have  much  in  common 
and  we  consider  them  in  parallel.   In  the  process  of  looking  for  a  better 
algorithm,  we  analyzed  and  compared  all  known  algorithms  for  finding  all 
circuits.   To  avoid  misunderstanding  we  will  define  our  terminology. 

A  graph  G  is  a  triple  <V,E,f>  where  V  is  a  finite  set  of  vertices, 
E  a  finite  set  of  edges,  and  f  is  a  function.   A  graph  is  either  directed 
(a  digraph)  when  f  :E  -*■  Vx  V,  or  is  undirected  when  f  :E  ->  { {u,v}  |u,v£V}  . 
Observe  that  this  definition  permits  the  graph  to  have  parallel  edges  (two 
edges  e.  and  e  are  parallel  if  e  ^   e  and  f(e  )  =  f(e  ))  and  self-loops 
(an  edge  e  is  a  self-loop  if  for  some  vertex  v,  f(e)  =  (v,v)  for  digraph 
or  f(e)  =  {v}  for  undirected  graph).   A  graph  G  is  s imple  if  it  does  not 
have  either  self-loops  or  parallel  edges.   An  edge  e  of  a  digraph  is  said 
to  be  incident  out  of  vertex  v  ,  and  incident  into  a  vertex  v  if 
f(e)  =  (v  ,v  ).   The  edge  e  is  incident  with  both  vertices  v,  and  vQ. 
The  in-degree  of  a  vertex  v  of  a  digraph  is  the  number  of  edges  incident 
into  v,  and  the  out -degree  of  v  is  the  number  of  edges  incident  out  of  v. 
An  edge  e  of  an  undirected  graph  is  incident  with  vertices  v  and  v  if 
f(e)  =  {v  ,v  }.   The  degree  of  a  vertex  v  of  an  undirected  graph  is  the 
number  of  edges  incident  with  this  vertex,  self -loops  being  counted  twice. 


Two  vertices  v  and  v  are  adjacent  if  there  is  some  edge  e  so 

that  f(e)  =  (v  ,v  )  or  (v  ,v  ,).   An  edge -sequence  from  vertex  v  to  v,  in  a 

digraph  is  an  alternating  sequence  of  vertices  and  edges,  v  eve  v  ... 

v  , e,  v,  ,  such  that  for  1  ^  i  S  k,  f(e.)  =  (v.  -,,v.).   This  edge-sequence 
k-1  kk;  i      i--L  i 

is  said  to  pass  through  the  vertices  v. ,  0  <  i  <  k.   An  edge-sequence  is 

simple  if  it  does  not  pass  through  a  vertex  more  than  once.   A  dipath  from 

v  to  v,  is  an  edge-sequence  from  v_  to  v,  where  all  vertices  are  distinct, 
Ok  Ok  ' 

except  possibly  that  v„  =  v  .   When  v  =  v,  ,  the  dipath  is  called  a  dicircuit. 
The  length  of  a  dipath  is  the  number  of  edges  in  it.   Since  the  vertices  can 
be  determined  knowing  the  edges  by  means  of  the  function  f,  vertices  may  be 
deleted  from  a  representation  of  an  edge-sequence.   Also,  since  we  will  be 
considering  an  edge-sequence  as  a  subgraph,  the  order  of  these  edges  is 
unimportant  and  an  edge-sequence  (and  hence  a  dipath  and  dicircuit)  is  simply 
a  set  of  these  edges  of  the  digraph. 

By  deleting  an  edge  e  from  a  graph  G  =  <  V,E,f  >  we  obtain  a  new 
graph  G'  =  <V,E',f  >  where  E'  =  E  -  {e}  and  f  is  a  restriction  of  f  to 
E'.   By  deleting  a  vertex  v  from  a  graph  G,  we  obtain  a  new  graph  G'  =  <  V  ,E',f>  • 
where  V'  =  V  -  {v}  and  E'  =  (eeE|f(e)  ^  (v,u)  or  (u,v)  for  ueV)  and 
f '  :E'  ->  V  x  V  is  a  restriction  of  f  to  E'.   By  ignoring  the  direction  of 
an  edge  e  we  mean  that  the  ordered  pair  f(e)  ==  (u,v)  be  considered  as  a  set 
{u,v}.   The  corresponding  definitions  for  an  undirected  graph  follow  if  the 
directions  of  the  edges  are  ignored. 

An  undirected  graph  is  connected  if  there  exists  a  path  between 
every  pair  of  vertices.   A  vertex  is  a  cut -vertex  if  it  lies  on  every  path 
between  a  pair  of  vertices.  ■  A  connected  undirected  graph  is  separable  if  it  has 
at- least  one  cut-vertex.   A  maximal,  connected  subgraph  of  a  graph  is  called  a 
component  of  the  graph.   A  maximal,  nons'eparable  subgraph  of  a  graph  is  called 
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a  block  of  the  graph.   A  directed  graph  is  strongly-connected  if  there 
exists  a  dipath  between  every  pair  of  vertices.   A  maximal,  strongly- 
connected  subgraph  of  a  digraph  is  called  a  fragment.   A  spanning  tree  of  a 
connected  undirected  graph  is  a  maximal  subgraph  which  has  no  circuits.   The 
unique  circuit  obtained  by  adding  a  non-tree  edge  to  the  spanning  tree  is 
called  a  fundamental  circuit  (with  respect  to  that  spanning  tree). 

Let  n  be  the  number  of  vertices  and  m  the  number  of  edges.  We 
shall  refer  to  an  edge-sequence  of  length  k  as  a  k-sequence.   Similarly,  we 
use  k-dipath  and  k-dicircuit. 


2.   Various  Methods  of  Circuit  Finding 

The  algorithms  for  finding  all  circuits  of  a  graph  fall  into  four 
broad  classes  depending  on  the  underlying  approach. 

1.  Circuit  vector  space  for  an  undirected  graph  (Maxwell 
and  Reed  (1965),  Welch  (1966),  Rao  and  Murti  (1969)) 

2.  Search  algorithms  (Char  (1968),  Chan  and  Chang  (1969), 
Roberts  and  Flore s  (1966),  Floyd  (1967),  Tier nan 
(1970),  Weinblatt  (1972)) 

3.  Powers  of  adjacency  matrix  (Ponstein  (1966),  Kamae 
(1967),  Yau  (1967),  Danielson  (1968)) 

h.      Edge-digraph  (Cartwright  and  Gleason  (1966)). 

The  details  of  these  algorithms  often  obscure  the  basic  ideas  behind  them. 

To  gain  insight,  and  to  be  able  to  point  out  specific  sources  of  inefficiencies, 

we  give  a  general  description  of  these  four  methods  of  circuit  finding.   Some 

new  results  are  also  given;  for  example,  it  is  shown  that  there  are  only  four 

undirected  graphs  whose  circuit  vector  space  consists  of  all  circuits  (ignoring 

the  null  element).   We  also  show  that  there  are  classes  of  graphs  for  which 

the  ratio  of  number  of  circuits  to  all  the  vectors  in  the  circuit  vector  space 
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goes  asymptotically  to  zero,  as  the  number  of  vertices  increases.   Unifying 
these  observations,  we  suggest  improvements.   A  new  class  of  undirected 
graphs  --  circuit-graphs  of  an  undirected  graph  —  are  defined  which  aid  the 
process  of  finding  circuits  in  undirected  graphs.   Some  properties  of  these 
graphs  are  explored. 

In  the  rest  of  the  paper,  by  the  term  "graph"  we  shall  refer 
simultaneously  to  undirected,  as  well  as  directed  graphs.   The  term  "circuits" 
will  also  be  used  likewise.   Should  there  be  a  possibility  of  confusion,  an 
explicit  mention  will  be  made  as  to  whether  we  are  discussing  the  directed 
or  undirected  graphs. 


2.1  Circuit  Vector  Space  of  an  Undirected  Graph. 

It  is  well  known  that  the  set  of  all  circuits  and  edge-disjoint 
unions  of  circuits  of  an  undirected  graph  form  a  vector  space  over  GF(2). 
The  vector  operation  is  the  "ring-sum"  operation  (see  Harary  (1969)  or  Deo 
(197^))-   A  non-null  vector  in  this  circuit  vector  space  is  called  a  circ. 
To  obtain  all  circuits,  one  finds  a  basis  for  the  circuit  vector  space 
(hereinafter  called  a  set  of  basic  circuits).   Let  p.   be  the  dimension  of  the 
space  (ju  is  also  called  the  nullity  of  the  graph).   All  the  2  -jLt-1  circs 
(excluding  the  basic  circuits  themselves,  and  the  null  element)  are  computed. 
Since  not  every  one  of  these  circs  is  a  circuit,  a  test  is  necessarily  made 
to  determine  if  the  circ  is  a  circuit. 

In  fact,  we  will  prove  that  there  are  only  four  graphs  whose 
circuit  vector  space  consists  of  just  circuits.   It  is  useful  to  have  the 
following  definition. 


-5- 

Definition  1:  A  graph  G  is  said  to  "be  reduced  if 

1.  G  is  simple 

2.  G  has  no  vertex  of  degree  0  or  1  and 

3.  if  there  is  any  vertex  v  in  G  of  degree  2,  then 
the  two  vertices  adjacent  to  v  are  adjacent. 

A  ring-sum  of  a  1-circuit  and  any  other  circuit  is  always  an  edge- 

'P\ 
disjoint  union  of  circuits.   Consider  the  set  P  of  all  (p   2-circuits  formed 

out  of  a  set  of  p  parallel  edges.      The  ring-sum  of  any  subset  S  of  these 

circuits  is  either  a  circuit  in  P  or  an  edge-disjoint  union  of  circuits  in 

P.   No  circuit  passes  through  a  vertex  of  degree  0  or  1.   If  a  circuit 

contains  an  edge  e  incident  with  a  vertex  v  of  degree  2,  it  must  also 

contain  the  other  edge  e  incident  with  v.   These  considerations  motivated  the 

above  definition. 

Observe  the  following  (see  Figure  l): 

1.  If  a  reduced  graph  is  not  K,,  the  complete  graph  on  three 
vertices,  then  the  nullity  \i   of  the  graph  is  at  least 
two. 

2.  There  is  no  reduced  graph  of  one  or  two  vertices. 
There  is  only  one  reduced  graph  of  three  vertices, 
viz.  ,  K~.   And  only  two  reduced  graphs  are  possible 
with  four  vertices  namely  K,  and  K.  -x.   In  all  these 
graphs  K~,  K,  and  K.  -x  every  circ  is  a  circuit. 

3.  A  reduced  graph  with  five  vertices  and  at  most  seven 
edges  must  have  an  edge-disjoint  union  of  circuits. 

LEMMA:  Let  G  be  a  reduced  graph  with  the  number  of  vertices,  n  ^  6.  If  G 
has  no  edge-disjoint  union  of  circuits,  then  the  number  of  edges  m  >  n + 3- 
The  converse  is  not  true. 
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K, 


K, 


K 


3,3 


Figure  1:   Graphs  all  of  whose  circs  are  circuits 


(a) 


00 


Figure  2:      Graphs  with  m   =  n+3  and  yet  having 
edge-disjoint  unions   of  circuits 
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PROOF:   Since  G  is  reduced,  and  has  no  edge-disjoint  union  of  circuits,  G 
must  be  connected.   No  vertex  in  G  can  be  of  degree  less  than  3;  for  otherwise 
G  is  either  not  reduced  or  has  an  edge-disjoint  union  of  circuits.   Therefore 
e^-=^=n  +  —  ^  n  +  3*  since  n  1  6. 

The  graphs  in  Figure  2  have  m  =  n+3  and  yet  have  edge-disjoint 
unions  of  circuits. 

THEOREM  1:   Only  the  four  reduced  graphs  K  ,  K.  ,  K.  -x  and  K    have  all  circs 
as  circuits. 

PROOF:  The  theorem  is  true  for  graphs  with  at  most  five  vertices  and  seven 

edges  by  the  observations  made  earlier.   It  is  true  for  graphs  with  n  vertices 

n  ^  6,  and  m  edges,  m  ^  n  +  2  by  the  lemma.   For  graphs  with  n  l  5,  and  m  1  n  +  3; 

the  theorem  will  be  proved  by  contradiction;  it  will  be  assumed  that  there  exists 

a  graph  G  which  has  n  5  5>  and  is  not  K   ,  and  yet  has  all  circs  as  circuits 

3  >  3 

and  draw  a  contradiction.   If  G  is  disconnected,  we  are  done.   Therefore,  G 
is  connected,  and  is  either  planar  or  nonplanar. 

Case  1:  G  planar,  n  ^  5,  m  S  n  +  3,   that  is  nullity  ju  =  m-n+1  is  at  least  four. 

There  are  at  least  /i  finite  regions,  and  an  infinite  region 
defined  on  a  plane  representation  of  the  graph  G   (see  Harary  (1969))- 
If  any  two  finite  regions  are  not  adjacent,  the  circuits  defining 
these  regions  are  edge-disjoint.   Likewise,  every  finite  region  must  also  be 
adjacent  to  the  infinite  region;  otherwise  the  circuit  defining  the  finite 
region  is  disjoint  with  the  circuit  defining  the  infinite  region.  The  dual 
G  is  therefore  a  complete  gr  ph  K    of  ji+1  vertices,  and  /i+l;-'5  which  is 


impossible.   Thus,  there  must  exist  a  pair  of  regions  which  are  not  adjacent, 
and  hence  an  edge-disjoint  union  of  circuits. 

Case  2:  G  is  nonplanar,  and  therefore  G  either  has  a  proper  subgraph 
homeomorphic  to  K  _,  or  has  a  subgraph  homeomorphic  to  K  . 

Case  2.1:   Let  G  have  a  subgraph  homeomorphic  to  K  .   An 
edge-disjoint  union  of  circuits  of  K  shown  in  Figure  3 
is  homeomorphic  to  a  subgraph  of  G. 

Case  2.2:  G  has  a  proper  subgraph  g  homeomorphic  to  K   . 

——————  3  )  J 

Since  g  is  a  proper  subgraph,  there  exists  an  additional  path 
between  two  vertices  u,  and  v.   If  u  and  v  correspond  to 
vertices  in  the  same  independent  set  of  vertices  of  K   , 
there  exists  an  edge-disjoint  union  of  circuits  (see  Figure  If  (a)). 
If  u  and  v  correspond  to  vertex  u'  in  one  independent  set  of 
vertices  of  K    and  to  v'  in  the  other  independent  set, 

3>  5 

respectively,  an  edge-disjoint  union  of  circuits  exists  (see 
Figure  1+  (b ) ) . 

Observe  that  the  Theorem  1  does  not  give  us  information  as  to  how 
many  edge-disjoint  unions  of  circuits  a  given  graph  can  have.   In  general,  a 
large  fraction  of  the  circs  are  edge-disjoint  unions  of  circuits.   The  ratio 
of  circuits  to  all  vectors  in  the  vector  space  for  numerous  classes  of  graphs 
tends  to  zero.   Three  such  classes  are  shown  in  Figure  5.   For  wheel  graphs, 
the  nullity  fi   is  n-1,  the  number  of  circuits  is  c  =  ju(jU-l)+l.   For  the  modified 
ladder  graphs,  ju  =—  (n+2)  and  c  =  ^-  +  -^  -  3.   For  complete  graphs 
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Figure  3:   An  Edge -disjoint  Union  of  Circuits 

in  a  Complete  Graph  of  Five  Vertices 


(a) 


00 


Figure  k:     Edge -disjoint  Unions  of  Circuits  in 
a  Nonplanar  Graph  Having  a  Proper 
Subgraph  Homeomorphic  to  K   . 
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(a)  A  Wheel  Graph 


4 4 4 


(b)  A  Modified  Ladder  Graph 


(c)  A  Complete  Graph 


Figure  5 
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2  n 

\i  =  ~  -   —^  +  .1,  and  c  =  -   Z  (g)  (s-l)!   In  all  the  three  cases  lim 

s  =3  n-»oo 

—  ->  0.   Following  these  assertions  we  make  the  following  conjecture: 

In  any  large  reduced  graph  the  ratio  of  circuits  to  the  circs  goes 
asymptotically  to  zero  with  the  number  of  vertices,  n. 

To  give  a  numerical  example  consider  K,»,  the  complete  graph  on  10 
vertices.   The  number  of  circuits  is  556  014  and  the  number  of  edge-disjoint 
unions  of  circuits  is  3k   359  192  353«   Thus  the  efficiency  of  an  algorithm 
for  finding  all  circuits  would  improve  dramatically  if  the  number  of  edge- 
disjoint  unions  of  circuits  it  computes  can  be  reduced.   An  attempt  to 
compute  only  a  subset  of  all  circs,  and  yet  find  all  circuits  was  made  by 
Welch  (1966).   He  chose  the  basis  as  a  set  of  fundamental  circuits,  and  tried 
to  give  an  ordering  of  the  fundamental  circuits,  namely  that  if  there  exist 

three  fundamental  circuits  C,  C.  and  C  ,  i  <  j,  and  C   H  C.  /  0,  and 

1   j      &  k    i 

Cv  DC.  ^  0,  but  C.  D  C.  =  0,  then  k  <  j.   Such  an  ordering,  if  existed, 
k    J  1    j 

would  have  reduced  the  number  of  edge-disjoint  unions  of  circuits  computed. 
Three  years  later  Gibbs  (1969)  gave  a  counterexample  and  proved  that  there 
are  graphs  for  which  such  an  ordering  does  not  exist.   The  corrected  Welch- 
Gibbs  algorithm  computes  all  ^   linear  combinations.   We  will  shortly  propose 
a  method  which  reduces  the  number  of  edge-disjoint  unions  of  circuits  computed 
in  the  process  of  finding  all  circuits. 

The  dicircuits  of  a  digraph  do  not  constitute  a  vector  space  under 
any  operation  (see  Williams  (1972)).  Thus  this  approach  cannot  be  used  with 
digraphs. 
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2.2  Search  Algorithms  ; 

There  are  two  kinds  of  search  algorithms.   Algorithms  due  to  Char 

(1968)  and  Chan  and  Chang  (1969)  produce  all  permutations  of  i  vertices  chosen 

out  of  n  vertices  of  a  graph  (l  ^  i  ^  n).   Any  graph  G  of  n  vertices  can  be 

considered  on  a  subgraph  of  the  complete  graph  K  .   A  circuit  of  G  is  also 

a  circuit  of  K  .   Thus  each  circuit  of  K  is  tested  to  see  if  it  contains  an 
n  n 

edge  not  present  in  G. 

The  second  kind  of  search  algorithms  construct  a  dipath  between 

vertices  v.  and  v.  and  test  to  see  if  there  is  an  edge  e  such  that  f(e)  =  (v., v.), 
1      J  J   1 

We  shall  describe  both  kinds  of  algorithms  in  detail  in  Section  3« 


2.3  Powers  of  Adjacency  Matrix 

Let  A  be  the  variable  adjacency  matrix  of  a  digraph.   Then  a  non- 
null  element  A  .  .,  p  1  1  is  the  set  of  all  p-sequences  from  vertex  v.  to 
ij  1 

vertex  v..   Since  no  dicircuit  can  be  of  length  greater  than  n,  we  need 
J 

compute  only  up  to  A  .   Observe  that  an  edge-sequence  need  neither  be  a  dipath 
nor  a  dicircuit.   The  crux  of  the  problem  is  to  avoid  such  (non-simple)  edge- 
sequences.   It  is  in  this  avoidance  that  different  algorithms  using  the  method 


-x- 
Def.  Fl.   A  variable  adjacency  matrix  A  of  a  digraph  D  =  <V,E,f>  is  a  nxn 

symbolic  matrix  in  which  the  (i,,i )  -element  A.  .  =  Z   e,  ,  where  e.  eE  such  that 

1J  all  ev  k        k 
k 

f(e  )  =  (v., v.).   If  there  is  no  such  en  ,  then  A.  .  =0.   The  powers  of  A,  p  ^  1, 

are  defined  as  follows:   A  =  A  and  for  p  >  1,  Ap  =  Ap   ■  A  =  A  .  aP"1 
employing  the  usual  symbolic  multiplication. 
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.V.  -r\ 

vary.   Kamae  (1967)  eliminates  dicircuits  and  flowers  of  length  p  from  A 
and  proceeds  to  the  next  power  of  A.   Others,  such  as  Danielson  (1968), 
Ponstein  (1966)  do  not  delete  the  circuits  and  flowers,  but  determine  whether 
or  not  a  given  edge-sequence  is  a  dicircuit. 


2.k     Edge  Digraph  of  a  Digraph 

Consider  a  digraph  £(D)  =  <J,Q,^>  derived  from  the  given  digraph 
D  =  <V,J,f>  as  follows:   The  vertex-set  of  g.(D)  is  the  edge-set  of  D,  and  the 
edge-set  Q,  of  £(D)  contains  an  edge  q,  ty(q)    =   (e  ,e  ),  e   and  e   in  J  iff 
for  some  v,  v  and  v  in  V,  f (e  )  =  (v  ,v)  and  f (e  )  =  (v,v  ).   Such  a  graph 
is  called  the  edge  digraph  of  D.   (See  Harary  (1969)) 

Now  consider  a  p-sequence  S  of  a  digraph  D'.   The  p-sequence  S 

yields  a  (p-l)- sequence  in  £  (D1),  a  (p-2)-sequence  in£(£(D')),  and  in  general, 

k         k-1 
a  (p-k) -sequence  in  £  (D')  =  (g  '  (D')),  for  k  ^  p.   Suppose  now  that  all 

1-dicircuits  are  reported  and  deleted  from  D',  resulting  in  a  subdigraph  D. 

Every  edge  in  £(D)  will  then  correspond  to  a  simple  2-sequence,  and  vice  versa. 

An  edge  e  of  £(D)  will  correspond  to  a  rooted  2-dicircuit   if  it  is  an  edge  of 

a  2-dicircuit  of  £(D).   Report  and  delete  all  2-dicircuits  of  £(d),  denoting 

the  resulting  digraph  as  £  (D).   (£.(D)  is  the  edge  digraph  of  £  ."  (d)  with 

all  j-dicircuits  deleted  from  it;  £. (D)  =   D. ) 


Def.  F2.   A  flower  of  a  digraph  is  a  subdigraph  which  is  a  union  of  a  dipath 
and  a  dicircuit  of  length  at  least  two  such  that  the  last  vertex  of  the  dipath 
is  in  the  dicircuit,  and  this  vertex  is  the  only  vertex  in  the  dicircuit  and 
dipath. 

Def.  F3.   A  rooted  k-dicircuit  is  a  dicircuit  of  length  k,  with  one  of  its 

k  vertices  specified  as  the  root.   Every  k-dicircuit  has  k  rooted  k-di circuits, 
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Assume  inductively  that  there  is  a  one-to-one  correspondence  between  - 
the  simple  p-sequences  of  D  and  the  edges  of  £  "  (D).   Report  and  delete  all 
its  p-dicircuits,  and  denote  the  resulting  digraph  as  t        (D).   Then  each  edge  in 
£P-1(D)  corresponds  to  a  p-dipath  of  D.   Let  £p(D)  =  £(£P~  (D)).   Now,  each  edge 
in  6P(D)  corresponds  to  two  adjacent  edges  of  £   (D),  or  what  is  the  same, 
corresponds  to  two  p-dipaths  P  and  P  of  D  having  (p-l)  edges  in  common. 
These  two  dipaths,  therefore,  form  a  simple  (p+l) -sequence  of  D.   Obviously,  if 
the  first  vertex  of  P  (or  P  )  is  the  same  as  the  last  vertex  of  P   (or  P..  ), 
their  union  gives  a  (p+l)-dicircuit  D-,;  and  D  yields  (p+l)  rooted  (p+l)-dicircuits, 
each  of  which  has  a  corresponding  edge  in  a  (p+l)-dicircuit  of  £/  (D) .   If  the 
edge  e  of  £  (D)  does  not  lie  in  any  (p+l)-dicircuit,  then  e  corresponds  to  a 
(p+l)-dipath  of  D. 

Thus,  given  a  digraph  D,  we  report  all  1-dicircuits,  delete  them  and 
proceed  to  take  the  edge  digraph  of  D;  report  and  delete  all  2-dicircuits  of 
the  edge  digraph  of  D,  and  so  on.   Observe  that  a  p-dicircuit  will  be  reported 
p  times  if  an  algorithm  is  implemented  using  the  method  as  described. 

Computationally  it  will  be  advantageous  to  perform  some  initial  simpli- 
fication on  the  given  graph  before  actually  proceeding  to  find  all  circuits. 
Parallel  edges,  self-loops,  vertices  of  degree  0  (in-degree  or  out-degree  0,  for 
digraphs)  vertices  of  degree  2  (for  undirected  graphs),  vertices  of  in-degree  and 
out-degree  1  can  all  be  taken  care  of  initially.   Then  if  the  given  graph  G  is 
undirected,  decompose  it  into  blocks.   If  G  is  a  digraph  decompose  it  into 
fragments.   The  circuit  finding  algorithms  can  now  be  applied  on  these  subgraphs. 
This  is  advantageous  because  such  editing  algorithms  have  a  time  bound  of 
0(max(m,n))  (see  Hopcroft  and  Tarjan  (1973)),  whereas  the  circuit  finding 
algorithm  is  essentially  an  exponential  algorithm  and  does  not  fall  in  any  of 
the  classes  defined  by  Karp  (1972)  and  Cook  (1971). 
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3.   Analyses  of  Algorithms 

Before  we  proceed  to  derive  time  bound  T.T  ' s  for  the  algorithms, 
let  us  keep  in  mind  the  caveat  that  as  far  as  possible  an  analysis  of  the 
'algorithm'  and  not  an  analysis  of  an  'implementation'  of  the  algorithm  is 
attempted.   However,  an  analysis  of  an  algorithm  without  any  regard  to  an 
implementation  is  not  always  possible.   So  much  so,  whenever  a  particular 
implementation  is  known  to  improve  the  speed  of  an  algorithm,  we  assume  that 
implementation  and  analyse.   We  assume  that  all  the  data  references  made  by  an 
algorithm  are  to  data  elements  stored  in  a  random  access  memory  and  that  once 
a  circuit  is  found  it  is  output  to  an  external  storage  medium  (and  also  stored 
simultaneously  in  the  random  access  memory  if  the  algorithm  requires  at  a 
later  stage  circuits  computed  earlier). 


3-1  Circuit  vector  Space  Algorithms 

Welch  (1966)-Gibbs  (1969)  algorithm:  The  Welch-Gibbs  algorithm  chooses 
a  set  of  fundamental  circuits  as  a  basis  for  the  circuit  vector  space.   The 
algorithm  as  given  by  Gibbs  (1969)  would  be  reporting  some  edge-disjoint  unions 
of  circuits  as  circuits  and  possibly  missing  some  circuits  if  a  set  of  non- 
fundamental  but  basic  circuits  were  used  as  a  starting  set. 

Let  B  =  fB,  ,  B  ,  . . .  ,  B  1  be  the  starting  set  of  fundamental  circuits. 

1    2  jU 

We  build  up  iteratively  two  sets  S  and  Q,  such  that  when  the  process  terminates 
S  contains  all  circuits  and  Q  contains  all  circs.   Initially  both  S  and  Q 
contain  fundamental  circuit  B,  only.   At  the  beginning  of  the  i   iteration 
2  ^  i  ^  ju,  S  and  Q  contain  all  circuits  and  circs  obtained  so  far.   During  the 
i^h  iteration  we  obtain  circuits  and  circs  by  taking  ring-sum  of  B.  with  each 
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element  M  of  set  Q,  and  then  add  them  to  sets  S  and  Q,  respectively.   If  M  and 
B.  do  not  intersect,  their  ring-stun  R  will  be  an  edge-disjoint  union  of  circuits, 
and  vector  R  is  included  in  a  set  Q.  which  contains  edge-disjoint  unions  of 

■4-Vi 

circuits  computed  during  the  i   iteration.   If  M  and  B.  do  intersect,  their 
ring-sum  R  could  be  a  circuit  or  an  edge-disjoint  union  of  circuits,  and  then 
vector  R  is  included  in  a  set  S,  which  contains  the  ring-sums  of  a  circ  M  and 
the  fundamental  circuit  B.  if  M  intersects  with  B. .   Thus,  circuits  computed 
during  the  i   iteration  are  elements  of  the  set  S, .   It  can  be  shown  that  a  ' 
vector  ¥  in  S  corresponds  to  a  circuit  iff  it  does  not  contain  or  is  contained 
in  another  vector  in  S  .   Vector  W  is  chosen  from  S,  ,    and  compared  with  every 

"U  ~Xj 

element  in  S  .   If  for  any  vector  R'  in  S  ,  W  <=  R'  then  R'  is  not  a  circuit, 
and  the  vector  R'  is  deleted  from  S  and  added  to  Q,  .   If,  on  the  other  hand, 
W  =3  R'  then  W  is  not  a  circuit,  and  the  vector  W  is  deleted  from  S  and  added 
to  G^_ .   Choose  a  new  vector  from  S,  as  W  and  repeat  the  process.   At  the  end 
of  this  search  for  containment,  S,  will  contain  all  and  only  circuits  computed 

th 

during  the  i°  iteration.   The  set  S,  and  the  fundamental  circuit  B.  are  added 

t  1 

to  both  S  and  Q,  set  Q,  is  also  added  to  Q,  and  we  proceed  to  the  next  iteration. 

Observe  that  the  set  of  circs  Q,  contains  exactly  2  -1  vectors  at  the 

th 
end  of  the  i   iteration.   Considering  a  worst  case  when  every   vector  in  Q,  will 

be  intersecting  with  the  fundamental  circuit-vector  B.  , ,  the  set  S,  of  all 

ring- sums  of  a  vector  M  in  Q  and  B.    such  that  M  and  Q  intersect,  will  contain 

at  most  21-!   vectors.   Let  EL  be  the  number  of  comparisons  made  in  the  testing 

for  the  minimality,  with  respect  to  containment  of  the  vectors  in  S, ,  during 

th 
the  i   iteration.   Then, 

N.  =  a.   |S,  I2   where   0  <  a.  <1 

i    l  '  t '  l 

If  S't  contained  3.  (2X-l)  vectors,  0  <  p.  <  1,  then 

N.  =a.  p2(2x-l)2 
1    1  Kl v    ' 
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Numbers  a.    and  p.  depend  on  the  structure  of  the  graph.   The  extreme  values 
are  not  likely  to  be  true  for  any  graph.   And  the  total  number  of  comparisons 

made  is   Z  a.    p2(2X-l)2  =  0(22^).   Thus,  the  time  bound 
i=2  1 
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All  the  2-1  vectors  computed  have  to  be  stored,  and  the  storage 
required  is  of  the  order  of  m  •  2   bits. 

Algorithm  of  Rao  and  Murti  (1969):   The  R  &  M  algorithm  differs  from 
W-G  algorithm  in  that  that  a  vector  which  is  a  ring-sum  of  some  basic  circuit- 
vectors  is  not  tested  for  its  minimality  with  respect  to  containment  of  another 
vector.   Instead  an  attempt  is  made  to  construct  a  spanning  tree  of  the  given 
graph  G  using  all  but  one  edges  present  in  a  circ.  We  will  succeed  in  this 
attempt  iff  the  circ  is  a  circuit.   Observe  that  it  is  not  necessary  to  have 
in  the  memory  all  the  circuits  and/or  edge-disjoint  unions  of  circuits  computed 
so  far,  and  that  every  circ  will  have  to  be  computed  as  a  ring- sum  of  some 
basic  circuits.   If  a  set  of  fundamental  circuits  is  chosen  as  the  basis,  and 
if  the  chords  are  labeled  from  1  to  \x,   then  we  need  only  store  the  characteristic 
submatrix  B,   of  the  fundamental  circuit  matrix  Bf  (see  Deo  (1970)-   The  matrix 
B,   is  a  binary  matrix  of  size  \x   x  (n-l). 

Let  R  be  the  ring-sum  of  k  fundamental  circuits,  and  let  it  contain 

I   edges.   Without  loss  of  generality,  R  =  {c,  ,  c  ,  ...,  c,  ,  b,,  b  ,  ...,  b  ,  ) 

where  c. 's  stand  for  the  chords  and  b. 's  for  tree  branches  with  respect  to  a 
11 

spanning  tree  T.   If  R  were  a  circuit,  a  spanning  tree  containing  any  (i-l) 
edges  of  R,  and  in  particular  a  spanning  tree  containing  all  the  tree  branches 
of  R  and  the  (k-l)  chords  c.  for  2  %   i  <   k,  must  exist.   On  the  other  hand,  if 
R  is  an  edge-disjoint  union  of  circuits,  it  properly  contains  a  circuit  which 
could  be  made  a  fundamental  circuit  with  respect  to  some  spanning  tree.   Thus, 
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the  chord  c  is  added  to  the  spanning  tree  T,  and  a  search  is  made  to  see  if   ; 
there  exists  a  tree  branch  of  T  not  in  R  that  could  he  deleted  from  T  to  break 
the  fundamental  circuit  containing  the  chord  c  .   If  there  is  no  such  branch 
b  available,  obviously  R  properly  contains  this  circuit,  and  R  is  not  a  circuit. 
Otherwise,  the  chord  c  is  introduced  into  T  and  b  is  deleted  from  T.   Call  the 
resulting  spanning  tree  as  T.   The  set  R  now  contains  one  less  a  chord  than  it 
did  prior  to  the  transformation.   This  process  is  repeated  for  all  chords  c. , 
2  ^  i  =  k,  if  necessary.   If  all  these  chords  could  be  made  as  tree  branches, 
then  R  is  a  circuit,  otherwise  it  is  not. 

All  2-1  vectors  are  computed  and  tested  to  see  if  a  vector  corresponds 
to  a  circuit.   A  vector  is  computed  as  a  ring-sum  of  some  basis  vectors.   The 
number  of  ring-sum  operations  performed  is  given  by 

k=2   k 

Given  a  vector  R,  a  ring-sum  of  k  basis  vectors,  testing  if  it  corresponds  to 
a  circuit,  in  the  worst  case,  needs  the  introduction  of  all  k-1  chords.   The 
total  number  of  chords  used  is,  again,  ju  •  2    -  \i  .   Thus  for  the  Rao  and  Murti 
algorithm  we  have  the  upper  bound 

TR&M  =  0^)  (3-1.2) 

It  is  easy  to  see  that  the  lower  bound  is  also 

Algorithm  of  Maxwell  and  Reed  (1965):  All  the  2^-1  circs  are  computed 
first;  then  each  circ  R  is  compared  with  every  other  circ  to  see  if  R  is  a 
circuit.   A  worst  case  analysis  would  immediately  give  an  upper  bound  of 
0(2  )  which  is  also  an  upper  bound  for  Welch-Gibbs  algorithm.   However,  direct 
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comparison  shows  that  for  any  graph  G,  Welch-Gibbs  algorithm  makes  lesser 
number  of  comparisons  for  containment  than  Maxwell  and  Reed  algorithm. 

Thus,  of  the  three  circuit  vector  space  algorithms  described  Rao 
and  Murti  algorithm  is  the  fastest,  and  requires  the  least  storage.   (However, 
see  Section  k>  ) 


3.2  Search  Algorithms 

These  algorithms  are,  in  general,  applicable  to  both  undirected,  as 
well  as  directed  graphs  unlike  the  circuit  vector  space  algorithms  which  are 
applicable  only  to  undirected  graphs. 

Char's  Algorithm  (1968):   The  algorithm  has  a  list  of  all  circuits, 
a  master  circuit  matrix,  of  a  complete  graph  K  .   If  all  circuits  of  a  graph 
G  having  n  vertices,  n  ^  p,  are  required,  each  circuit  of  K   of  length  at 
most  n  is  tested  to  see  if  it  contains  an  edge  not  in  G.   If  all  the  edges  in 
the  circuit  considered  are  also  edges  of  G,  the  circuit  is  reported.   Assuming 
that  the  circuits  of  a  complete  graph  K.,  j  ^  n  are  listed  lexicographically 

J 

earlier  than  the  circuits  of  the  complete  graph  K .  -, ,   the  time  bound  is 

Tc  =0(|   Z   (*)-(s-l)!)  (3.2.1) 

s  =3 

Storage  is  required  for  all  the  circuits  of  the  complete  graph,  K  . 

Assuming  each  circuit  is  stored  in  one  word,  the  total  requirement  is 

n 

Z  (  )  (s-l)!  words.   Clearly,  storage-wise  Rao  and  Murti  algorithm  is  vastly 
s=3 

superior  to  Char's  algorithm.   Since  no  closed  form  expression  is  available 
for  Eq. (3.2.I),  it  is  difficult  to  compare  Char's  algorithm  with  Rao  and  Murti 
algorithm.   Assuming  that  the  sum  on  the  right-hand  side  of  Eq. (3. 2.1)  is 
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proportional  to  nn ,   Rao  and  Murti  algorithm  would  be  superior  to  Char's  when 


i.e.,  jti  +  %2jU  g  nfc^  n. 

That  is,  Rao  and  Murti  algorithm  is  superior  to  Char's,  when  the  nullity  \i   is 

linear  in  the  number  of  vertices,  or  equivalently,  if  the  number  of  edges  m  is 

linear  in  n.  .. 

Algorithm  of  Chan  and  Chang  (1969):   The  Chan  and  Chang  algorithm  uses 

a  modified  variable  adjacency  matrix.   Actually,  it  generates  all  permutations 

and  tests  if  there  exists  a  directed  edge  between  two  consecutive  vertices  in 

the  permutation  and  between  a  chosen  vertex  and  the  first  vertex.   Given  an 

initial  sequence  1  2  3  •  •  •  n  1>  ^e   permutation  P  is  pivoted  at  the  r-th 

element  p  in  P,  1  <  r  ^  n.   All  integers  p. ,  i  ^  r,  are  successively  incremented 

by  1.   The  number  n  when  incremented  results  in  r,  by  definition.   Once  a  new 

permutation  is  thus  obtained  every  subsequence  1  2  3  • •  sis  tested  to  see  if 

a  dicircuit  exists  with  1_2_3_. -_s_l  as  its  vertex- sequence.   We  generate  all 

permutations  containing  1  and  then  delete  1  and  keep  2  as  the  first  vertex  in 

the  next  set  of  permutations.   In  general,  once  all  permutations  with  k  as  the 

first  vertex  are  generated,  delete  k  and  take  k+1  as  the  first  vertex  and 

proceed.   Terminate  if  k=n.   The  time  bound  for  this  algorithm  is 

n-1 
Tc&c  =0(  Z  i!)  (3.2.2) 

i  =  3 

Algorithm  of  Roberts  and  Flores  (1966),  Floyd  (1967)  and  Tiernan 
(1970):  The  algorithm  which  was  originally  given  by  Roberts  and  Flores  (1966), 
for  finding  all  Hamiltonian  dicircuits  of  a  digraph,  in  actuality  generated  all 
dicircuits.   In  Floyd  (1967),  it  appears  as  a  nondeterministic  algorithm. 
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Tiernan  (1970)  gives  a  deterministic  version  of  the  algorithm,  and  proves  its 
finiteness  and  correctness.   The  algorithm  as  given  in  Roberts  and  Flores 
(1966)  and  Tiernan  (1970 )  uses  the  successor  listing  .   Every  vertex 
v  of  the  digraph  has  a  row  in  a  matrix,  and  all  the  end-vertices  of 
edges  incident  out  of  the  vertex  v  are  listed  as  the  successors  of  v.   But 
improvement  in  the  speed  occurs  if  the  adjacency  matrix  of  the  digraph  is 
used. 

The  algorithm  constructs  a  dipath  P  and  tests  if  there  is  a  directed 
edge  from  the  last  vertex  to  the  first  vertex  of  the  dipath  P.   If  there  is 
such  a  directed  edge  a  dicircuit  is  found,  and  is  duly  reported.   The  construc- 
tion of  a  dipath  is  as  follows:   Let  p  be  the  first  vertex  of  the  dipath  being 
constructed.   If  there  is  a  directed  edge  (p  ,  p  )  such  that  p  >  p  ,  add  the 
edge  to  the  dipath,  and  -che  length  of  the  dipath  increases  by  1.   Suppose  that 
a  k-dipath  P  is  constructed,  so  far.   Let  its  vertex- sequence  be 
P  =  P-,-pG-..-p,   -p,  .   An  attempt  is  now  made  to  extend  the  dipath  P.   This 
extension  is  possible  only  if  there  is  an  edge  (p  ,  v)  such  that  v  is  not  a 
vertex  in  P,  because  otherwise  a  flower  blooms!   The  extension  is  made  only  if 
v  is  not  'forbidden'  and  v  >  p  .   This  device  ensures  that  every  dicircuit  is 
reported  once  only  by  defining  a  one-to-one  and  onto  mapping  of  dicircuits  to 
rooted  dicircuits,  rooted  at  its  lowest  numbered  vertex.   If  no  directed  edge 
(pv>  v)  satisfying  the  above  conditions  exists,  a  check  is  made  to  see  if  a 
dicircuit  can  be  found  as  described  earlier.   If  a  circuit  is  confirmed,  it  is 
reported.   In  any  case,  the  last  vertex  of  the  dipath  p  is  deleted.   Call  the 
subdipath  p  -p  -. . . -p,    as  P,  and  the  above  extension  process  is  now  applicable 
to  the  new  dipath  P.   The  last  vertex  p,  is  not  available  for  extension, 
because  otherwise  a  dipath  already  reported  would  result.   Thus,  the  vertex 
pv  is  made  'forbidden'  from  the  vertex  p   .   This  does  not  necessarily  mean 

K.  A  —  -L 
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that  p  is  forbidden  from  other  vertices,  or  that  p,  will  be  forbidden  forever 
k  -& 

from  p  ., .   Once  the  last  vertex  (which  is  pn  ,  )  of  the  new  dipath  Pis 
k-1  k-J. 

similarly  deleted,  the  vertex  p,  is  made  available  from  p,   . 

The  algorithm  starts  with  p  =  1.   A  new  first  vertex  p,  «-  p,  +  1  is 
used  after  the  current  first  vertex  p,  is  deleted.   The  algorithm  terminates 
after  detecting  that  the  vertex  p,  =  n. 

Tiernan  (1970 )  claims  that  his  algorithm  is  the ' "theoretically  most 
efficient"  search  algorithm.   The  claim  has  to  be  disputed  for  two  reasons: 

1.  Tiernan ' s  algorithm  does  consider  a  dicircuit  more 
than  once  during  the  search. 

2.  Even  a  search  algorithm  which  considers  each 
dicircuit  once  and  only  once  cannot  be  said  to 
be  theoretically  most  efficient. 

To  show  that  Tiernan' s  algorithm  does  consider  a  dicircuit  more  than  once,  we 

exhibit  the  graph  in  Figure  6.   The  dicircuit  1-2-3-1  will  be  found  earlier, 


Figure  6:  A  Graph  for  the  Tiernan' s  Algorithm 

but  will  not  be  reported  until  the  dipath  1-2-3-1+  is  explored.   Also,  in  testing 
for  a  directed  edge  from  vertex  k-   to  vertex  1,  a  directed  edge  (k,   2)  is  found 
giving  the  dicircuit  2-3-I+-2;  yet,  this  dicircuit  will  not  be  reported  until 
the  first  vertex  of  the  dipath  becomes  2.   The  above  observations  are 
implementation  independent,,  i.e.,  whether  a  successor  listing  is  used  or  a 
binary  adjacency  matrix  is  used  these  dicircuits  are  found  once  before  they  are 
found  again  and  reported. 
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Secondly,  in  trying  to  consider  a  di circuit  once  only,  Tiernan's 
algorithm  considers  each  directed  edge  a  number  of  times.   Minimizing  the 
number  of  examinations  of  each  individual  edge  also  improves  the  efficiency. 
Thus,  a  more  efficient  algorithm  might  be  possible  by  minimizing  the  number 
of  considerations  of  each  dicircuit  in  conjunction  with  a  minimization  of 
the  number  of  examinations  of  each  edge. 

To  estimate  a  lower  bound  for  the  RFFT  algorithm,  observe  that  for 
each  dipath  considered,  a  test  will  be  made  to  see  if  a  directed  edge  exists 
from  the  last  vertex  of  the  dipath  to  the  first  vertex,  and  the  list  containing 
the  forbidden  vertices  from  the  last-but-one  vertex  of  the  dipath  is  cleared. 
The  list  is  implemented  as  a  binary,  linear  array;  n  operations  are  required 
to  clear  the  list.   Now,  the  algorithm  considers  all  dipaths  and  only  dipaths 
with  the  property  (which  for  short  we  shall  call  property  jt),  that  the  first 
vertex  is  the  least  numbered  vertex  in  the  dipath.   Thus,  the  time  bound  is 

0(n  •  (number  of  dipaths  with  property  n))  =  T___m      (3«2.3) 

Since  for  each  k-dicircuit  there  are  at  least  (k-l)  dipaths  with  the  property 
tt,  the  algorithm  is  not  linear  with  respect  to  the  number  of  dicircuit s,  as 
claimed  by  Tiernan  (1970). 

Weinblatt's  Algorithm  (1972):  Weinblatt  suggests  an  algorithm  similar 
to  the  RFFT  algorithm.   A  dipath  P  is  built  up,  and  if  there  exists  an  edge 
from  the  last  vertex  of  the  dipath  P  to  another  vertex  in  P,  a  dicircuit 
thus  formed  is  immediately  reported,  in  contrast  to  the  RFFT  algorithm.   Also, 
the  Weinblatt's  algorithm  considers  each  directed  edge  once  only,  and  makes  use 
of  parts  of  dicircuits  reported  earlier  to  construct  new  dicircuits.   For 
example,  if  a  dipath  P  has  its  last  vertex  on  a  dicircuit  D, ,  a  test  is  made 
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to  see  if  any  other  vertex  of  P  is  on  the  dicircu.it  D, .   Or,  if  the  dicireuit 
D  had  a  vertex  in  common  with  another  dicireuit  D  which  in  turn  had  a  vertex  on 
the  dipath  P,  these  two  dicircuits  D  and  D  are  used  to  give  a  new  dicireuit. 
In  general,  D  ,  D  ,  . . . ,  D,  may  be  used  similarly  to  give  a  new  dicireuit. 
However,  this  necessitates  a  test  to  see  if  a  part  of  the  new  dicireuit  as 
given  by  D  ,  D  ,  . . .  ,  D  is  also  not  given  by  another  set  of  dicircuits  D.  , 

D.  ,  ... ,  D.  . 

2         J 

The  algorithm  is  recursive  and  giving  a  non-trivial  bound  on  the 
depth  of  recursion  is  a  difficult  combinatorial  problem,  and  the  analysis  of 
the  algorithm  is  thereby  complicated.   We  will  not  therefore  attempt  to 
analyze  the  algorithm.   However,  it  appears  to  be  faster  than  the  RFFT 
algorithm. 


5.3  Algorithms  Using  Powers  of  Adjacency  Matrix 

The  algorithms  will  be  described  first,  and  then  all  of  them  will 
be  directly  compared  at  the  end  of  this  subsection.   (Note  that  the  Def.  F2 
for  A  ,  powers  of  variable  adjacency  matrix  assumes  that  no  entries  in  A 
are  modified  as  will  be  done  in  the  following  algorithms.   Hence,  in  these 
algorithms,  A    •  A  £  A  •  A   ,  in  general.  ) 

Yau's  Algorithm  (1967):  A  nonzero  entry  A  ij  gives  all  the  p- 
sequences  from  vertex  v.  to  vertex  v.  as  a  sum  of  products,  each  product 
giving  a  p-sequence.   If  a  p-sequence  is  non-simple,  then  it  must  contain  a 
k-dicircuit,  for  k  <  p.   Yau's  algorithm  computes  A  and  tests  each  p-sequence 
in  a  non-zero  entry  of  A  to  see  if  it  contains  any  k-dicircuit,  k  <  p,  reported 
earlier.   If  it  does,  that  p-sequence  is  deleted  from  the  sum.   Obviously,  if  a 
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p-sequence  does  not  contain  any  k-dicircuit,  k  <  p,  it  must  "be  simple.   If 
the  p-sequence  is  from  a  non-diagonal  entry,  it  is  a  dipath;  otherwise  it 
is  a  dicircuit  and  hence  is  reported.   The  diagonal  entries  in  A  are  set 
to  zero  once  all  the  p-dicircuits,  if  any,  are  reported,  and  p  is  incremented 
by  1  and  A  *-  A    •  A,  and  the  process  is  repeated.   The  algorithm  terminates 
if  either  p  =  n  or  if  for  some  p,  A  is  an  all- zero  matrix.   This  termination 
clause  is  common  to  all  algorithms  described  in  this  subsection,  and  we  shall 
not  repeat  it. 

Ponstein's  Algorithm  (1966) :   If  A  contains  all  simple  p-dipaths 
and  only  those,  then  we  can  construct  A    satisfying  the  same  property  as 
follows:  Compute  A  #A  and  A^  -A.   Include  in  A . .   only  those  (p+l)-sequences 
which  are  present  in  (A*  A  ).  .  and  (A  •  A).  ..   Assume  that  a  (p+l)-sequence  R 

-'-J  -^  J 

in  (A  •  A  ) .  .  is  non-simple.   Since  A  contained  only  simple  edge-sequences,  R 

is  non-simple  because  of  an  edge  (v.,  v ),  where  v.  was  already  an  interior 

vertex  on  the  dipath  from  v  to  v..   Observe  that  in  the  (p+l) -sequence  R, 

k     J 

the  only  vertex  of  out-degree  2  is  v. .   No  (p+l)-sequence  in  (A  •  A).  .  can  be 
an  edge-sequence  with  v.  having  an  out-degree  of  2,  because  this  can  happen 
only  if  A?,  were  non-null.   Thus  R  will  not  appear  in  A  •  A.   It  is  easy  to  see 
that  a  simple  p-sequence  S  will  appear  in  both  A  •  A  and  A  •  A.   The  diagonal 
entries  give  rooted  (p+l)-dicircuits.   Report  and  delete  these  dicircuits. 
Obviously,  A    now  satisfies  the  stated  property. 

Danielson's  Algorithm  (1968):  A  major  source  of  inefficiency  in 
Yau's  algorithm  is  in  the  testing  to  see  if  an  edge-sequence  is  simple. 
Danielson's  algorithm  employs  a  much  simpler  test. 

All  edges  incident  into  a  vertex  v  are  labeled  with  v.   Using  these 
labels  construct  the  variable  adjacency  matrix  A.   Construct  another  matrix  A  , 
obtained  from  A  by  replacing  every  nonzero  entry  by  the  null  string  X.      Any 
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string  s  concatenated  (multiplied)  with  \   gives  the  same  string  s.   Call  as    ; 
A2  the  product  A  •  A  ,  and  AP  =  A  •  AP~  ,  p  £  3.   Then  an  edge-sequence  appearing 
in  AP  also  gives  corresponding  vertex- sequence  but  with  its  end-vertices 
missing.   Assuming  that  AP~  contains  only  (p-l)-dipaths,  a  p-sequence  from 

vertex  v  to  vertex  v.  is  non-simple  iff  v.  appears  in  the  p-sequence.   All 

1  J  1 

non-simple  edge- sequences  are  deleted  from  A ,  and  all  p-dicircuits  are 
reported  and  deleted.   Obviously,  A  now  contains  only  p-dipaths.   Observe 
that  because  of  this  peculiar  edge-labeling,  only  premultiplying  of  A    with 
A  is  allowed. 

The  above  three  algorithms  report  the  rooted  dicircuits.   If  a 
dicircuit  need  be  reported  only  once,  we  may  report  the  unique  rooted 
dicircuit  of  every  dicircuit  rooted  at  the  smallest  numbered  vertex  in  the 
dicircuit. 

Kamae's  Algorithm  (1967):   This  algorithm  uses  the  binary  adjacency 

matrix,  X.   Powers  of  this  matrix  are  computed  in  real  field;  thus,  a  non-zero 

entry  X.T  gives  the  number  of  edge-sequences  from  vertex  v.  to  vertex  v.. 
ij  oh  1  j 

Assume,  as  before,  that  X    gives  the  number  of  (p-l)-dipaths  only.   Then 
to  compute  Xp  giving  the  number  of  p-dipaths  and  to  report  all  p-dicircuits 
during  such  computation,  it  will  be  necessary  to  find  the  number  of  non-simple 
edge -sequences.   If  the  matrix  X  were  to  be  obtained  by  postmultiplying  X 
by  X,  all  non-simple  edge-sequences  will  be  flowers  (see  Def.  F2. ).   Therefore, 
it  will  be  necessary  to  compute  the  number  of  p-flowers  from  vertex  v.  to 
vertex  v.  and  then  subtract  it  from  (XP~  •  X).  ..   This  number  is  the  sum  of 
the  number  of  flowers  with  k-dicircuits  (containing  the  vertex  v.)  and  with 

J 

(p-k)-dipaths  from  vertex  v.  to  vertex  v.,   2  ^  k  ^  p-1.   Every  k-dicircuit 
containing  v.  but  not  v.  is  considered,  and  the  number  of  (p-k)-dipaths  from 
Vi  to  Yi  is  comPuted.   Now,  the  length  of  the  dicircuit  k  takes  values  from  2 
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to  (p-l),  and  this  computation  is  performed  for  these  values  of  k  and  the 
sum  of  the  number  of  flowers  found  for  each  k  gives  the  total  number  of 
flowers  from  vertex  v.  to  vertex  v..   Construct  a  flower  matrix  F,  such  that 

J 

F. .  gives  the  number  of  flowers  from  vertex  v.  to  vertex  v..   Then 

-'-J  d 

P      "D-1 

X*  «-  Xy       •  X  -  F 
Now,  X?.  gives  the  number  of  simple  p-sequences  from  v.  to  v..   A  non-zero 

•*-d  ■!-        d 

diagonal  entry  X..  gives  the  number  of  p-dicircuits  containing  the  vertex  v. . 
Once  these  dicircuits  are  found  and  reported,  set  the  diagonal  entries  to 
zero,  increment  p  by  1  and  repeat  the  procedure. 

To  find  the  p-dicircuits,  compute  a  matrix  Z  defined  as  follows: 

z. .  =  x1?:1  •  x. . 
ij      ji       ij 

The  above  equation  self-explains  that  Z.  .  will  be  nonzero  iff  there  is  a 
p-dicircuit  containing  the  directed  edge  (v.,  v.).   For  details  of  tracing 
the  dicircuits  using  matrix  Z,  and  for  the  details  of  computing  the  flower 
matrix  F,  the  reader  is  referred  to  Kamae  (1967). 

The  algorithms  of  Danielson  (1968),  Kamae  (1967),  Ponstein  (1966), 
and  Yau  (1967)  will  now  be  directly  compared  with  each  other.   As  was 
mentioned  earlier  in  Section  2,   these  algorithms  differ  only  in  identifying 
the  non-simple  sequences,  and  it  is  this  identification  that  we  compare. 

Yau's  algorithm  identifies  a  non-simple  p-edge-sequences  by  testing 
it  to  see  if  it  contains  any  k-dicircuit,  k  <  p,  reported  earlier.   Thus, 

at  most  N  ,  tests  are  performed,  where  N  ,  is  the  number  of  dicircuits  of 

p-1  ^       '       p-1 

length  at  most  (p-l).   Ponstein  finds  a  non-simple  p-edge-sequence  of 

(A    •  A).  .  by  testing  to  see  if  the  sequence  is  absent  in  (A  •  A   ).  .,  and 

-*-d  -*-d 

this  would  require  at  most  s  tests,  where  s  is  the  number  of  p-edge-sequences 
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in  (A*  AP   ).  ..   All  flowers  from  the  vertex  v.  to  vertex  v.  are  computed 
by  Kamae's  algorithm  and  this  requires  considering  all  k-dicircuits,  k  <  p, 
containing  v.  hut  not  v..   For  each  k-dicircuit  considered,  the  edges  of 

J-  J 

incident  out  of  the  vertices  of  the  dicircuit  are  deleted,  and  dipaths  in  the 
resulting  digraph  connecting  the  vertex  v.  with  vertex  v.  are  to  he  counted. 

J 

This  would,  in  turn,  require  that  the  flowers  of  the  resulting  digraph  be 
enumerated,  and  so  on.   Thus,  Yau's  algorithm  is  faster  than  Kamae's. 

Danielson's  algorithm  determines  whether  or  not  a  given  p-sequence 
from  vertex  v.  to  vertex  v.  is  simple  by  testing  if  the  vertex- sequence 
contains  the  vertex  v..   Only  one  operation  is  required  for  each  sequence. 
Since  only  premultiplying  was  allowed,  any  non-simple  edge-sequence  is  a 
flower.   The  time  bound  is  therefore 

T  =  0 (Number  of  all  dipaths  +  Number  of  all  flowers).     (3.3.I) 

For  each  sequence  Ponstein's  algorithm  makes  as  many  tests  as  there  are 
terms  in  (A  *A ). .,  and  Yau's  algorithm  makes,  at  most,  as  many  tests  as  there 
are  previously  reported  dicircuits.   Danielson's  algorithm  is  therefore 
superior  to  those  of  Yau,  Ponstein  or  Kamae. 

Comparing  with  RFFT  algorithm,  the  set  of  dipaths  with  property  n 
is  a  subset  of  all  dipaths,  and 

TD  -  TEFFT  *  (5. 3- 2 J 


3- 2. k     Edge  Digraph  Algorithm 


The  concept  of  edge  digraph  is  used  in  an  algorithm  by  Cartwright 
and  Gleason  (1966).   The  method  was  described  in  detail  in  Section  2.k   and  we 
will  not  repeat  it  here. 
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The  union  of  two  simple  p-dipaths  P  and  P  will  be  a  (p+l)- 

dicircuit  if  the  first  vertex  of  P..  is  the  same  as  the  last  vertex  of  P 

1  2 

and  there  is  a  common  (p-l)-siibdipath  in  P,  and  P  .   Thus,  to  find  a  (p+l)- 

dicircuit,  a  pair  of  p-dipaths  P,  and  P  having  a  (p-l)-dipath  in  common, 

and  also  having  a  last  vertex  of  P  as  the  first  vertex  of  P  must  be  found. 

Thus,  for  each  (p+l),  N  •  N  p-sequence  pairs  are  considered.   So  a  lower 

bound  is 

n-1 

°(  Z.    #  ■  TC&G  <5^l) 

P=l 

where  N  is  the  number  of  simple  p- sequences. 

The  number  of  simple  p-sequences  is  not  less  than  the  number  of 
dipaths  with  the  property  tt  of  the  RFFT  algorithm.   Moreover  Eq.  (3.V.I)  is 
quadratic  in  N  . 

The  RFFT  algorithm  for  dicircuits  is,  therefore,  the  fastest  of  all 
algorithms  analyzed,  and  Rao  and  Murti  algorithm  is  the  best  of  the  three 
circuit  vector  space  algorithms  analyzed.   Weinblatt's  algorithm  appears  to 
be  faster  than  Tiernan's  algorithm  for  the  reason  that  parts  of  some  circuits 
computed  earlier,  are  used  to  produce  new  circuits. 


k-      Circuit-graph 

As  was  already  pointed  out  in  Section  3. 1  the  main  source  of 
inefficiency  in  the  algorithms  using  the  circuit  vector  space  approach  is  in 
the  wasted  computation  of  edge-disjoint  unions  of  circuits.   It  would  be  ideal 
if  we  computed  only  the  circuits.   Toward  this  end  we  define  a  new  graph  — 
a  circuit -graph  of  a  graph  --  and  demonstrate  its  use  in  the  computation  of 
all  circuits  of  a  graph. 
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Let  B  =  {B,,  B  ,  . . .  ,  B  )  be  a  set  of  "basic  circuits.   Observe 
the  following: 

1.  Any  circuit  C  of  a  graph  G  can  be  expressed  as  a 
ring-sum  of  b  basic  circuits,  1  ^  b  ^  \i, 

C  =  B.   ©  B.   ©  ...  ©  B.  . 
Xl    X2  Xb 

Consider  this  set  of  basic  circuits,  S  =  {B.  ,  B.  , 

1    2 

. ..,  B.  ).   It  is  not  possible  to  partition 

S  into  two  nonempty  blocks  S,  and  S  such  that  the 
ring- sum  C-,  of  all  circuits  in  S,  is  disjoint  with  the 
ring-sum  C  of  all  circuits  in  S  . 

2.  If  two  circs  C,  and  C  do  not  intersect,  i.e.  if 

C-,  D  C   =  0  then  their  ring- sum  C  ©  C  cannot  be  a 
circuit.   However,  if  C,  and  C  do  intersect,  C  ©  C 
may  or  may  not  be  a  circuit. 

Definition  2:  A  circuit-graph  of  a  graph  G  =  <V,E,f>  with  respect  to  a  set 

B  of  basic  circuits  of  G  is  a  graph  ^(G,B)  =  <V  ,E  ,f  >.   The  vertex  set  V  =  B, 

and  for  an  edge  eeE  ,  f„(e)  =  {B. ,  B.)  iff  B.  ^  B.  and  B.  D  B.  ^  0  for  1  gi,j  gjj 

An  example  of  a  circuit-graph  of  GT ,  the  modified  ladder  graph  is 

given  in  Figure  7-   The  set  of  fundamental  circuits  with  respect  to  the 

spanning  tree  shown  in  bold  lines  is  chosen  as  B.   For  a  different  set  of 

basic  circuits  we  would,  in  general,  get  a  different  circuit-graph.   But  the 

number  of  vertices  of  a  circuit-graph  of  G  is  always  jii,  the  nullity  of  G. 

Observe  that  if  the  graph  G  has  p  blocks,  then  every  circuit-graph  (with 

respect  to  any  set  of  basic  circuits)  of  G  will  have  p  components.   Consider 

a  subset  S  =  {B.  ,  B.  ,  . . .  ,  B.  }  of  vertices  of  a  circuit-graph  6(G,B)  of  G. 
Xl  X2  Xb 

If  B  =  B.   ©  B.   ©  ...  ©  B.   is  a  circuit  of  G,  then  there  must  be  a  path  in 
12  b 

C(G,B)  between  any  two  vertices  B.  ,  and  B.   in  S.   Thus,  with  any  circuit  C 

Xk       %  . 
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B, 


;a)   g. 


B„       B. 


(b)  A  Circuit-Graph  Gn   of  GT 

C         Li 


Figure  7. 


o 
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f  G  we  can  associate  a  connected  subgraph  of  £(G,B)  with  vertices 


B    B   ,  ... ,  B.  .   Only  those  subsets  of  vertices  of  g(G^B)  which  are 
the  vertex-sets  of  connected,  subgraphs  of  £>(g,B)  are  of  interest.   Any 


representative  of  these  connected  subgraphs  (with  the  same  vertex-set)  may 
be  taken.   A  convenient  subgraph  is  the  connected  induced  subgraph.    This 
correspondence  between  circuits  of  G  and  connected  induced  subgraphs  of  £(G,B) 
is  one-to-one  but  not  necessarily  onto,  because  of  observation  2. 

To  compute  all  circuits  of  a  graph  G,  then  the  procedure  is  as 
follows : 

Procedure  CCUBE:   [Compute  all  Circuits  using  a  Circuit-graph]: 

1.  Take  a  set  B  of  basic  circuits  of  G. 

2.  Construct  the  circuit-graph  c!(G,B). 

J.   Find  the  vertex-sets,  V's  of  all  connected, 
induced  subgraphs  of  (^(G,B). 

Ij..   Compute  the  ring-sum  R  of  the  basic  circuits 

corresponding  to  the  vertices  in  the  vertex-set  V. 

5.  Test  if  R  is  a  circuit.   If  yes,  report  the 
circuit. 

6.  Take  a  new  vertex-set  V,  and  go  to  Step  Ij..   If 
no  new  vertex- set  V  exists,  stop. 

Observe  that  a  test  in  Step  5  to  decide  if  S  is  a  circuit  is 
necessary. 

Let  us  apply  this  procedure  to  GT  (Figure  7)-   The  circuit-graph 
Gc  produced  in  Step  2  is  given  in  Figure  7(b).   (The  basis  B  was  the  set  of 
fundamental  circuits  with  respect  to  the  spanning  tree  shown  in  bold  lines.  ) 


Def-  F1+:  An  induced  subgraph  g  =  <VT,E',f '>  of  a  graph  G  =  <V,E,f>  is  a 
maximal  subgraph  of  G  with  the  vertex  set  V . 
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(a)   GT  with  a  Hamiltonian  Path  as  a  Spanning  Tree 
h 


(b )  Circuit-graph  6(GL,B ' ) 

(B'  is  the  set  of  fundamental  circuits  with 
respect  to  the  spanning  tree  shown  in  (a)) 


Figure  8. 
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There  are  exactly  —  +  7^  -  3  connected  induced  subgraphs,  where  \±   is  the 

number  of  vertices  of  G^  (nullity  of  GT).   That  is  every  connected  induced 

subgraph  of  Gp  corresponds  to  a  unique  circuit  of  G,,  and  there  is  no  wasted 

computation.   However,  for  a  different  basis  this  property  may  not  hold. 

In  fact,  if  a  Hamiltonian  path  is  chosen  as  the  spanning  tree  to  produce  a 

set  B'  of  fundamental  circuits  thenIS  (GL,B' )  is  a  complete  graph  of  ju  vertices 

(see  Figure  8).   This  brings  us  to  an  important  question: 

Q:  What  basis  B  of  the  circuit  vector  space  of  the 

graph  G  be  chosen  so  that£(G,B)  has  a  minimum  number 
of  connected  induced  subgraphs? 

No  precise  answer  has  been  found.   Approximate  answers  to  this  question  are: 

Choose  that  basis  B  for  which 

Al:  the  sum  of  lengths  of  the  basic  circuits  is  a  minimum 

A2:  the  number  of  edges  ofC(G,B)  is  a  minimum. 
The  algorithm  of  Paton  (1969)  attempts  to  produce  a  set  of  fundamental  circuits 
satisfying  Al  provided  the  vertices  are  labeled  'properly',  and  no  criterion 
is  known  to  label  these  vertices  'properly'.   (See  also  Prabhaker  (1972).) 
However,  it  may  be  noted  that  Paton' s  algorithm  never  takes  a  Hamiltonian 
path  (if  it  exists)  of  a  graph  G,  as  the  spanning  tree  to  produce  the  set  of 
fundamental  circuits,  unless  of  course,  G  is  a  circuit.   That  is,  if  we  use 
this  algorithm  to  find  a  basis  B,  the  circuit-graph i6(G,B)  may  not  be  complete. 
But  are  there  graphs  for  which  every  circuit-graph  is  complete?  We  have  the 
following  theorem  which  is  a  consequence  of  Theorem  1. 

THEOREM  2:   If  a  reduced  graph  G  is  not  one  of  the  four  graphs  K_,  K.  ,  K.  -x, 

J         k         k- 

or  K    then  G  has  a  circuit -graph  which  is  not  complete. 


* 

We  are  ignoring  the  trivial  connected  induced  subgraphs  which  are  the 
individual  vertices  of  the  graph. 
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A  brute  force  algorithm  for  finding  the  sets  of  vertices  of 
connected  induced  subgraphs  of  a  graph  can  be  easily  given  incorporating 
a  connected  components  algorithm  (say  of  Hopcroft  and  Tarjan  (1973))  which 
will  have  a  time  bound  of  0(ju  •  2  ).   For  each  set  of  vertices  of  a  connected 
induced  subgraph  of  (f  (G,B),  only  at  most  \±   ring-sum  operations  are  required 
(Step  k) •      To  test  if  a  circ  is  a  circuit  again  requires  at  most  ju  operations 
(using  a  test  similar  to  that  made  by  Rao  and  Murti  algorithm).   Now,  assuming 
that  an  algorithm  for  finding  the  vertex-sets  of  all  connected  induced 
subgraphs  of  a  graph  has  a  time  bound  of  0(n),  (where  N  is  the  number  of 
connected  induced  subgraphs  of C(G,B))  the  algorithm  for  all  circuits 
of  G  based  on  CCUBE  will  have  a  time  bound  of  0(ju  •  N).   Because  of 
Theorem  2,  we  can  always  choose  a  circuit-graph  which  is  not  complete, 
and  hence  N  <  2  -  ju  -  1.   Then  our  algorithm  is  faster  than  that 
of  Rao  and  Murti  algorithm.   However,  no  algorithm  for  finding  a  basis  of 
circuit  vector  space  of  a  graph  G  and  answering  the  question  Q,  is  known  to 
the  authors.   Thus  the  success  of  an  algorithm  based  on  CCUBE  depends  on  our 
being  able  to  give  efficient  algorithms  for  finding  the  vertex-sets  of 
connected  induced  subgraphs  of  a  graph,  and  for  finding  a  'good1  basis  for  the 
circuit  vector  space  of  a  graph. 


5.   Conclusions 

Of  all  the  analyzed  algorithms  for  finding  all  dicircuits  of  a 
digraph,  the  RFFT  search  algorithm  is  the  fastest,  and  requires  least  storage. 
Weinblatt  (1972)  algorithm  is  not  analyzed  because  of  its  extreme  complexity 
due  to  the  recursion.   It  appears  that  Weinblatt' s  algorithm  will  be  faster 
than  RFFT  algorithm.   The  time  bound  for  RFFT  algorithm  is 

TDTTir,m  =  0(n  •  (number  of  dipaths  with  the  property  it)) 

Kb  r  1 
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The  algorithm  of  Rao  and  Murti  (1969)  is  the  fastest  of  the 
three  circuit  vector  space  algorithms  known.  The  time  "bound  for  this 
algorithm  is 

T    =  0 (jLt  -  2M)        (  u   is  nullity  of  the  graph  G) 

And  if  ju  is  proportional  to  n  ,   Char's  (1968)  master  circuit  matrix 
algorithm  is  faster. 

Using  a  circuit-graph  as  a  guide  in  finding  all  circuits,  the 
number  of  edge-disjoint  unions  of  circuits  considered  can  be  reduced.   The 
efficiency  of  an  algorithm  for  all  circuits  based  on  a  circuit-graph  depends 
on  the  efficiency  of  an  algorithm  for  finding  all  connected  induced  subgraphs 
of  a  graph.   We  have  proved  that  there  are  only  four  graphs,  viz. ,  K,,  K.  , 
K,  -x  and  K,  ,,  having  all  circs  as  circuits.   Thus  for  any  other  graph  G, 
a  circuit-graph  of  G  can  be  constructed  which  is  not  complete.   The  number 
of  connected  induced  subgraphs  will  be  less  than  2  .   When  a  brute-force 
algorithm  is  given  for  connected  induced  subgraphs,  the  algorithm  for  all 
circuits  based  on  circuit-graph  has  the  same  time  bound  as  that  of  Rao  and 
Marti. 
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